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Abstract 

J 

A  sonobuoy  field  is  a  random  two  dimensional  array  if  the  signals 
i  from  the  sensors  are  coherently  processed.  The  sensor  position 

coordinates  must  be  estimated  to  a  fraction  of  a  wavelength  in  order  to 
coherently  process  the  signals,  which  is  a  separate  processing  task  since 
the  field  has  a  changing  random  pattern  due  to  drifting  of  the  buoys. 

This  paper  reviews  the  statistical  properties  of  the  maximum  likelihood 
bearing  estimator  computed  from  coherent  summation  of  output  signals  from 
randomly  deployed  sensors.  It  is  shown  that  the  approximate  mean-square 
bearing  error  for  a  small  random  array  favorably  compares  with  that  for  a 
square  array  of  equally  spaced  sensors  with  about  one-half  the  aperture. 
The  random  array  is  not  subject  to  the  aliasing  (grating  lobes)  problem 
of  an  equally  spaced  array,  and  a  great  deal  of  array  gain  can  be 
achieved  from  a  large  array  of  randomly  placed  sensors.  Serious 
consideration  should  be  given  to  tracking  the  locations  of  drifting 
sonobuoys  so  as  to  coherently  process  them  as  an  array. 
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Bearing  Estimation  using  Random  Arrays 


Melvin  J.  Hlnich 

Introduction 

A  sonobuoy  field  is  a  random  two  dimensional  array  if  the  signals 
from  the  sensors  are  coherently  processed  (Thorn,  Booth,  and  Lockwood 
[ 1 ] ) .  The  sensor  position  coordinates  must  be  estimated  to  a  fraction  of 
a  wavelength  in  order  to  coherently  process  the  signals,  which  is  a 
separate  processing  task  since  the  field  has  a  changing  random  pattern 
due  to  drifting  of  the  buoys.  It  is  technically  feasible  to  track  the 
pattern,  but  the  processing  gain  from  such  an  effort  must  be  perceived  as 
sufficiently  great  to  justify  such  an  effort.  This  paper  reviews  the 
statistical  properties  of  the  maximum  likelihood  bearing  estimator 
computed  from  coherent  summation  of  output  signals  from  randomly  deployed 
sensors.  Stochastic  properties  of  the  array  response  function  are 
discussed  by  Thorn,  et.  al.  and  Steinberg  [2]. 

The  signal  received  at  a  sensor  from  a  distant  acoustic  source  in  a 
horizontally  homogeneous  ocean  is  a  plane  wave  plus  ambient  noise  (Clay 
and  Medwin,  (3]).  Using  complex  variable  notation,  a  single  frequency 
plane  wave  at  point  (xj,x2)  and  time  t  is 

p(t,x)  *  A  exp[i2ir  fQ(t-(xi  cos60  +  x2  sin90)/c)] 

where  0O  is  the  horizontal  direction  of  propagation  with  respect  to  the 
x-axis,  c  is  the  phase  velocity,  f0  is  the  wave's  frequency  (the 
wavelength  is  XD  *  c/f0)  and  A  is  the  complex  amplitude.  This  amplitude 
will  be  assumed  to  be  a  constant  for  the  expositon  of  coherent 


processing,  but  actual  signals  will  have  some  amplitude  and  phase 
modulation  which  reduces  the  signal  coherence  across  the  array.  A 
reduction  of  coherence  results  in  an  increase  in  the  mean  square  error  of 
bearing  computed  from  any  of  the  coherent  processing  methods  discussed 
next. 

1.  Coherent  Processing 

The  most  widely  used  coherent  processing  technique  is  delay-and-sum 

beamforming.  Suppose  that  the  signal  received  by  the  kth  sensor  in  a 

horizontal  planar  array  is  sCt.x^)  »  p(t,x)t)  +  e(t,xjc)  where  e(t,x)c)  is  a 

Gaussian  noise  field  at  time  t  and  point  (xnc»X2k.).  The  usual  approach 

for  beamforming  is  to  compute  the  beam  output  signal  B(t,0^)  *= 

M 

£fc-l8(t+Tk£)  for  a  set  of  delays  TkJl  =  (xikcosSjj,  +  x2ksl-neJl  )/c  for 
sensors  k»l,...,M  and  beam  (look)  angles  £=2,...,N.  The  energy  in  beam 
over  a  sampling  window  of  length  T  is  |B(0£>|2  where  B(0^)  is  the 
filtered  beam  signal  in  a  band  about  fQ  of  bandwidth  1/T.  If  the  plane 
wave  signal  is  broadband,  then  the  filter  passband  is  adjusted  to  the 
bandwidth  of  the  signal. 

Levin  [A]  shows  that  the  maximum  likelihood  estimator  of  the 

A 

propagation  angle  0O  is  the  0O  in  the  interval  0  <  0  <  2ir  that  maximizes 
|B(0)|2.  Define 

A  2  M  A  £  M 

°1  “  M”1  l  (xlk  “  xi)2,  02  -  M-1  '  -  X2)2  (1.1) 

k-1  k-i 

_  ,  M  _  M  M  _  _ 

where  xj  -  M"lE k-lxlk  and  X2  -  M"lj:k«lx2k*  Set  Ek-l^lk-*) (x2k~y)  *  0 
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A2  *2 

by  a  proper  rotation  of  the  coordinate  system.  Then  If  pMoj  and  pM02  are 
large,  where  p  denotes  the  power  signal-to-noise  ratio  In  a  narrow  band 
about  fQ,  then  the  mean-square  error  of  0O  is  given  by  the  approximation 
(Hinich  [5]) 

A  A  -2  A-2 
E(0O  ~  e0)2  ®  (A0/2x )2(2pM)~l(a i  sin20o  +  02  cos20o)  rad%.  (1.2) 

This  approximation  also  holds  for  a  discrete  set  of  beam  angles  and  an 
estimator  0^  that  maximizes  |B(0$)|2,  provided  that  the  error  |0£+i-0£| 
is  small  for  all  l .  In  many  applications,  however,  the  bias  due  to  the 
discreteness  of  the  0£  is  larger  than  the  square-root  of  the  right  hand 
side  of  (1.2).  It  is  important  in  these  applications  to  interpolate 

A 

between  0jq  and  its  adjacent  angles. 

A  simple  way  to  interpolate  is  to  use  the  frequency-wavenumber 
equivalence  with  beamforming  (Hinich  [6]).  Suppose  that  s(t,xjc)  is 
passed  through  the  same  narrow  band  filter  as  was  used  to  produce  B(0). 
For  example,  the  fast  Fourier  transform  (FFT)  can  be  used  to  filter 
discretely  sampled  signals  to  yield  the  filtered  output 

.  N_1 

s(*k)  *  s(nx  ,jyl)exp(-i2iimn/N) ,  (1.3) 

where  T  »  Nr  is  the  sampling  period,  r  is  the  sampling  interval,  and  m  is 
the  integer  part  of  f0T.  Then 
M 

B(8*  )  ■  l  s(xjc)expli(K  ijjxifc  +  «2ix2k)l  (1.4) 

k-1 

where  <  -  (2x  /Ao)cos0£  and  <2i  *  (2ir /Ao)sin0j[, .  The  FFT  can  be  used  to 

compute  (1.4)  for  as  fine  a  grid  as  desired  by  imbedding  {xjj}  in  a  square 
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lattice  of  equally  spaced  points,  and  adding  zeros  to  expand  the  two 
dimensional  array. 

From  now  on  it  will  be  assumed  that  the  bearing  interpolation  is 
sufficiently  refined  and  pM  is  sufficiently  large  so  that  expression 
(1.2)  holds  for  0O.  This  expression  will  now  be  used  to  analyze  the 
ex-ante  root  mean-square  error  of  80,  rmse  60,  when  the  sensor  locations 
are  randomly  determined. 

2.  Gaussian  Random  Arrays 

Suppose  that  the  x^  and  X2k  are  M  independent  realizations  from  two 

independent  Gaussian  random  variables  Xj  and  X2.  With  no  loss  of 

generality  center  the  origin  of  the  coordinate  system  so  that 
2  2  2 

X!  -  N(0,o  1)  and  X2  ~  N(0,O2)>  where  on  is  the  variance  of  Xn  for  n-1,2. 

To  further  simplify  the  analysis,  let  Oj  *  02  “  o. 

M  _  M  _ 

It  then  follows  that  a_2£k*i(xik  -  xj)2  and  o-2£k=i(x2k  -  X2>^  are 

2 

independent  chi-square  x  variates  with  M-l  degrees  of  freedom.  Denote 
2  M_1  2 

0T  *  Pr(x  >  T2) ,  e.g.  0.99  **  Pr(x  >  5.2)  using  a  standard  chi-square 
M—  1  1 5 

distribution  table.  By  independence,  the  probability  that 
M  M 

Ek-l(xlk  “  xl)2  >  (oT)2  ,  Sk*l(x2k  “  x2)2  >  (°T)2 

2 

is  Bt»  and  thus 

a  _2  a  2  2 

Pr{M~koi  sin20  +  o2  cos20 )  <  (oT)-2}  >  0T  .  (2.1) 

Applying  (2.1)  to  (1.2)  and  ignoring  the  approximation  error. 


_ _ _  ii  in  ir "  '  ...  .  _  I  -fill n I 
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mse  0O<  (X0/8irp  oT)2  rad2  (2.2) 

2 

with  probability  greater  than  67.  For  example  if  £5?  =  0.99,  M  -  16, 
a  »  6X0,  and  p  *  1/10  (-10  db) ,  then  T2  -  5.2  and  from  (2.2), 
rmse  60  <  0.53°  with  probability  greater  than  0.98. 

A 

It  is  useful  to  compare  this  result  with  the  mean-square  error  of  60 
for  a  square  lattice  array  whose  M=N2  sensors  are  at  the  points 
{md,nd:m,n=l , . . . ,N}  .  Setting  L  »  (N-l)d  for  the  length  of  the  square's 
sides, 

M  _  M 

l  (xik  -  xi)2  =  l  (x2k  -  X£)2  »  ML2/12.  (2.3) 

k-1  k-1 

A 

Thus  from  (1.2)  the  mse  0O  for  the  random  array  is  less  than  that  for  the 

2 

square  lattice  array  (with  probability  greater  than  if 
1/2 

L  <  (12/M)  To.  (2.4) 

For  example  if  M*16  and  By  *  0.99,  then  the  bound  is  L  <  1.98  o.  If  M  * 

36  and  Bt  *  0.99,  T2  -  18.5  and  L  <  2.48  a. 

If  the  square  is  centered  on  the  origin  rather  than  at 

( (tffl)/d,(N+l)/d) ,  the  probability  that  a  sensor  is  in  the  square  is 
2 

f  Pr( | Y |  <  L/2o)l  where  Y  -  N(0,1).  Thus  if  L  ■  4o ,  this  probability  is 
0.91.  Since  (0.99)2  ■  0.98,  it  follows  from  the  examples  above  that 
the  small  Gaussian  random  array  has  approximately  the  same  bearing 
accuracy  as  a  square  lattice  array  of  about  one-half  the  aperture. 
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This  comparison  with  the  lattice  array  has  been  made  very 
conservative  by  the  selection  of  S?  =  0.99.  The  rmse  0O  for  a  small 
random  array  will  often  be  as  small  as  a  lattice  array  of  equal 
aperture. 

3.  Conclusion 

It  Is  shown  that  the  approximate  mean-square  bearing  error  for  a 
random  array  favorably  compares  with  that  for  a  square  array  of  equally 
spaced  sensors  with  about  one-half  the  aperture.  The  random  array  Is  not 
subject  to  the  aliasing  (grating  lobes)  problem  of  an  equally  spaced 
array,  and  a  great  deal  of  array  gain  can  be  achieved  from  a  large  array 
of  randomly  placed  sensors.  Serious  consideration  should  be  given  to 
tracking  the  locations  of  drifting  sonobuoys  so  as  to  coherently  process 
then  as  an  array. 
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